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Spin chains have long been considered as candidates for quantum channels to facilitate quantum 
communication. We consider the transfer of a single excitation along a spin-1/2 chain governed by 
Heisenberg-type interactions. We build on the work of Balachandran and Gong [J], and show that 
by applying optimal control to an external parabolic magnetic field, one can drastically increase 
the propagation rate by two orders of magnitude. In particular, we show that the theoretical 
maximum propagation rate can be reached, where the propagation of the excitation takes the form 
of a dispersed wave. We conclude that optimal control is not only a useful tool for experimental 
application, but also for theoretical enquiry into the physical limits and dynamics of many-body 
quantum systems. 



> 
in 

m 

^' 

o 
o 



X 



I. INTRODUCTION 

Quantum computers promise to allow efficient simula- 
tion of large dynamic and complex systems and deliver 
performance advantages over tlieir classical counterparts. 
One of the central considerations for the construction of 
a quantum computer is an infrastructure that can rapidly 
and robustly transport qubit states between sites where 
qubit operations can be performed. The components for 
this infrastructure may be though of as quantum chan- 
nels for quantum information transfer. One of the tech- 
nologies under investigation to constitute such a channel 
is the one-dimensional spin-chain [l|-[lo|. which consists 
of a string of particles coupled via their spin degrees of 
freedom, each acting as an effective two-level quantum 
system. As is customary in quantum information pro- 
cessing, proper engineering of the control parameters of 
the system is essential to achieve the high fidelity nec- 
essary for robust quantum computation. This can be 
obtained, for instance, by employing a numerical opti- 
misation method which, for the specific settings of the 
problem, seeks the optimal control pulses that allow one 
to implement the desired operation [lT| - [2l| . In this pa- 
per, we apply such a method, known in the literature 
as the Krotov method [23426| . to the case of quantum 
state transfer along a one-dimensional spin chain. The 
specific system we use was introduced by Balachandran 
and Gong [l| , but here we show that by designing the ex- 
ternal driving parameters with optimal control methods, 
one can obtain a significant increase in fidelity, even over 
short time scales [27| . 

These high-fidelity, high-speed transmissions exhibit 
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interesting characteristics. If one ignores the effects tak- 
ing place near the boundaries, the evolution of the exci- 
tation is that of a dispersed wave, moving with almost 
constant velocity along the chain. This velocity is in- 
dependent of the chain length, and furthermore has an 
upper bound, indicating the presence of a fundamental 
limit on the rate of transmission. Through a closer anal- 
ysis, we show that this limit can be directly related to 
the theoretical maximum speed of the state transfer al- 
lowable by the laws of quantum mechanics [28l - l33j . 

Producing time-optimal gates has already been ex- 
plored in the literature [3J-|39| where the authors con- 
sidered geodesies on the Bloch sphere for systems with a 
low number of dimensions. Unfortunately, extension of 
these methods to many-body systems (such as the system 
we consider here) are prohibitively difficult. Conversely, 
the numerical optimisation methods that we employ have 
little difficulty in finding sets of optimal solutions, even 
at this limit. In effect, we demonstrate that through ap- 
plication of optimal control, we can not only transmit 
the excitation with a high fidelity, but also at the fastest 
possible speed. One can even reverse the problem, im- 
plying that optimal control can be used to probe such 
fundamental dynamical limits on many-body quantum 
systems. Such tools will be invaluable as the ambition 
of quantum science leads it towards investigations of sys- 
tems of greater complexity which are less tractable ana- 
lytically. 

The paper is arranged as follows. In Section |TT1 we de- 
scribe the system used for information transfer in more 
detail and the precise scheme which we will use for prop- 
agating quantum information in the system. Section IIIII 
discusses the application of optimal control to the trans- 
fer scheme, and shows that optimal control can effect sig- 
nificant gains in the transfer speed. We then discuss the 
fundamental limit of these improvements in Section IIVI 



and show that optimal control in fact allows us to reach 
this limit, thus allowing us to transfer the spin state in 
the fastest possible time allowable by the laws of quan- 
tum mechanics. Finally, we present the conclusions in 
Section El 
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II. SPIN CHAINS AS QUANTUM CHANNELS 

A. Overview 

Using spin chains as quantum channels for communi- 
cation between two parties was first proposed by Bose 
in 2003 @ and later developed in a series of papers (we 
refer the reader to Ref. [H for a review). The idea is rel- 
atively simple: Alice (the sender) has a quantum state 
she wants to relay to Bob (the receiver). Between them 
is a onc-dimcnsional chain of N spin- 1/2 particles which 
are coupled via nearest-neighbour interactions. Alice has 
access to the first spin in this chain, and can prepare its 
spin state as she chooses. Bob has access to the final site 
(the terms 'spin' and 'site' will be used interchangeably), 
whose state he can read out. Following we apply an 
external parabolic magnetic field, which Alice can con- 
trol. The procedure for sending quantum information 
along the chain is as follows. 

1. The spin chain is prepared in its ground state with 
respect to the external magnetic field. 

2. Alice prepares the initial spin state to be the state 
she wishes to transfer. 

3. By manipulating the magnetic field, Alice controls 
the propagation of the spin along the chain, which 
takes place due to the coupling between the spin 
degrees of freedom. 

4. After some prescribed time when the state has been 
transferred to the final site. Bob reads out the state 
of this site. 



B. The Hilbert space and Hamiltonian 

The model we consider is sketched in Fig. [TJ It is 
composed of a one-dimensional spin- 1/2 chain with N 
sites, where distances are measured by the variable x 
(although this may not be a physical distance). We will 
consider uniform Heisenberg nearest-neighbour couplings 
characterised by the same coupling strength J, and the 
presence of a parabolic external magnetic field in the z- 
direction, normal to the direction x. Consequently, the 
field will act on the nth site as 



B„{t) ^ C{t){xn - d{t)y 



(1) 



where d{t) is the position of the field minimum along x at 
time t, and C{t) is a measure of the global field strength. 



FIG. 1. (Colour online) The one-dimensional spin chain used 
for information transfer. The (blue) filled circles represent 
sites along the chain, with the applied magnetic field depicted. 
The effective couplings are indicated operating between the 
sites. 



The Hamiltonian then takes the form 

N-l N 



= E • '^"+1 + E Bnit)a'n : (2) 



n=l 



n=l 



where n labels the spin sites, with n = 1 and n = N 
referring to the first and last spins, respectively, and 
CT„ = (cr^, (tK, cr^) are the Pauli spin operators for the nth 
spin. For convenience, all system parameters are scaled 
to make them dimensionless, and the coupling strength 
is set to J = 1. 

The dynamics are governed by the interplay between 
the nearest-neighbour interactions and the interaction 
of each site with the external parabolic magnetic field. 
When sites arc far from the field minimum, the local 
field strength dominates over the nearest-neighbour inter- 
actions, effectively 'switching off' the coupling between 
sites. For sites near the minimum where the field is 
weak, the nearest-neighbour coupling dominates, and the 
neighbouring sites interact with each other. These two 
processes control the propagation of spin states along the 
chain. 



C. Communicating quantum information 

We identify the computational basis for our system 
with the quantised states of each spin, such that |0) = \].) 
(spin down with respect to z) and |1) = |t) (spin up). 
Assume that Alice prepares the chain in the initial state 
^'(0)), with the first spin site in the state |1}, and all 
other sites to in their ground state |0). We can write this 
state as 

|«'(i = 0)) = = |1) ® |0) ® |0) ® • • • ® |0) , (3) 

with the first spin site in the state jl), and all other sites 
to in their ground state |0). The states \(pn) are defined 
as 



N 



I \Smn) , n = 1, 



,,iV, 



(4) 



where Smn is the Kronecker delta. Alice's goal is to ma- 
nipulate the magnetic field parameters C (t) and d{t) such 
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FIG. 2. (Colour online) The transfer begins with the state 
l^'(O)) and with the potential minimum centred at a; = 0. 
(a) The excitation is localised at the first site, (b) The field 
minimum moves along the chain during the evolution, (c) 
The spin state has been completely transferred to the final 
site in the chain. 



that at the final time T the final state obeys 



1 



(5) 



The protocol for transferring the state is based on that 
described in Ref. [l| , which we outline in Figure [5] The 
transfer begins with the state |^'(0)) and with the poten- 
tial minimum centred at x = 0. At first, the interaction 
between the first two sites dominates over the interaction 
with the locally weak magnetic field, and so the sites in- 
teract and the spin state migrates from the first site to 
the second. As the field minimum moves along the a:-axis, 
nearest-neighbour interactions are effectively switched on 
for pairs of spins closest to the minimum, and switched 
off for spin pairs that are distant. By correctly moving 
the field minimum and adjusting the field strength, the 
spin state is able to traverse the chain. The condition in 
Eq. ([S]) means that we do not preserve the phase of the 
initial state; to achieve this (as also discussed in Ref.[l|) 
one can use dual-rail encoding [40l - |42j . whereby one en- 
codes the qubit in the entanglement phase of a pair of 
spin chains. In what follows, we shall only consider the 
phase-insensitive transfer of a single excitation. 

Since each site of the chain has two internal spin states, 
the size of the Hilbert space V. scales exponentially with 
the number of sites, so that for N{> 1) sites, dim'H = 2^. 

However, since lH{t),J2n=i '^n] = the state of the sys- 
tem only evolves within the subspace U C H, spanned 
by the A'' basis states |0„) [Ij. The reduced size of the 
effective Hilbert space is particularly beneficial when one 
wants to numerically simulate the evolution efficiently. 
We do this by solving the associated Schrodingcr cqua- 



i-\^it)) = Hit)\^{t)) 



(6) 



where H is the matrix form of the Hamiltonian that acts 
only on the subspace U 



H{t)=Ho+Hiit), 



(7) 



with 



Ho = -2 J + J 



(\ 1 
1 I 
1 



yo 



0\ 



1 

1 V 



(8) 



and 



Hr(t) = diag (/o(i), ■ • • , /jv-i W) , (9) 



where /^(t) = C (t) {xn — d{t)) (note that we have 
rescaled the energy so that spins pointing down do not 
contribute to the total energy). The Schrodingcr equa- 
tion is integrated numerically using the Crank-Nicolson 
method [H]. 

This scheme was first considered by Balachandran and 
Gong who showed that by choosing d{t) = st and 
C{t) = k, where s and k arc constant, one is able to 
adiabatically transfer the state across the entire chain 
with relatively good fidelities. However, the transfer rates 
here are very slow, with transfer times typically on the 
order of lO'' J for fidelities greater than 99%. 

In many proposed implementations of quantum com- 
puters, it is likely that transport processes will take up a 
significant amount of the total operating time. It there- 
fore seems clear that one should seek to minimise the 
time required for these processes. However, according 
to quantum mechanics there is some fundamental limit 
which restricts the speed at which we can communicate 
with our spin chain, referred to in the literature as the 
'quantum speed limit' (QSL) [2i,[2i,[31 33,44]. The goal 
is to come as close as possible to this limit, effectively 
communicating at the highest possible speed allowable 
by quantum mechanics. We shall see in the next section 
that optimal control can help us in this endeavour. 



III. OPTIMAL DYNAMICS 

We can state our problem in the following way: we 
start with an initial state, and want to control the system 
to produce the desired final state. In our case, the initial 
state is \(pi), and we want to achieve the final state \(Pn) 
(up to a global phase). We can control the evolution of 
the system using the external magnetic field, in particular 
the time-dependent controls d{t) and C{t). (Although in 
principle we could also control the inter-spin coupling J, 
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this is much more difficult to achieve experimcntaUy.) 
Optimal control theory provides us with a set of tools to 
search for the optimal way to control the system, often 
referred to as the set of optimal controls. 

Here, we implement an optimal control algorithm most 
commonly known as the Krotov method. In outline, the 
method works as follows. 

1. We solve the Schrodinger equation from ([6]) to find 
|5'(T)), where T is the total evolution time. 

2. We define the co-state |x(T)) = \^{T)){^{T)\tpN)- 
This state is propagated backwards to the initial 
time. 

3. The initial state is then propagated forward again 
through time, but at each time step we calculate 
the matrix elements 



(10) 



for the two controls ui{t) = d[t) and U2[t) = C{t). 
The matrix elements are then used to update the 
control functions, which are then used to propagate 
'^{t) to the next time step. 

4. We can then calculate the fidelity of the transport 

FEE|(vI/(^)|^^.)|^ (11) 

which tells us how close we were to achieving our 
goal. (Note that we will often refer not to the fi- 
delity, but to the infidelity / = \ — F.) If we achieve 
fidelity F = 1 (up to a given threshold), we stop 
the optimisation, otherwise we begin again at step 

m 

There are several aspects in implementing the alg orithm 
which are described in more detail in Ref. [22|. Fig- 
ure [3] shows the non-adiabatic transfer of a spin exci- 
tation across a chain of = 101 spins without apply- 
ing optimal control. One sees that during propagation, 
much of the spin excitation has been left behind. One 
way to correct this would be to lower the field strength: 
this will allow neighbouring sites to interact for longer, 
so that more of the excitation can be transmitted. How- 
ever, this causes the excitation to spread out, which can 
be seen in Fig. U) In comparison with Fig. [31 we see that 
although we have not left as much of the excitation be- 
hind, we have spread it over more sites. After applying 
optimal control (300 iterations of the update procedure), 
we arrive at the evolution shown in Fig. [SJ Here we see 
that we no longer leave excitation behind at the initial 
spin sites, and although we spread out the excitation dur- 
ing transport, we successfully recover the highly localised 
final state, giving a final fidelity F that differs from unity 
by < 10~^. The pulses required to achieve this result are 
shown in Fig. [HI Typical features of these pulses are large 
modulations at the boundaries, necessary for 'accelerat- 
ing' (it will be useful here to imagine an excitation wave) 
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FIG. 3. (Color online) Excitation probability plotted against 
X for a spin chain with = 101 sites, at times (a) 0, (b) 100 
and (c) 200, in units of Here, d{t) = 0.5£ and C{t) = 1. 
The final fidelity is only around 15%. 
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FIG. 4. (Color online) Excitation probability plotted against 
X for a spin chain with A*' = 101 sites, at times (a) 0, (b) 100 
and (c) 200, in units of Here, d{t) = 0.5t and C{t) = 0.1. 
The final fidelity here is around 5%. 
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FIG. 5. (Color online) Excitation probability plotted against 
X for a spin chain with TV = 101 sites, at times (a) 0, (b) 
100 and (c) 200, in units of J"^ Here, d{t) and c{t) were 
optimally controlled. The final fidelity is > 99%. 
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t = r/4, (c) t = r/2, (d) t = 3r/4, and (e) f = T, where 
r = IgsL = 56.50J"^ Both d(t) and C{t) were found after 
100,000 iterations of the optimal control algorithm. 



FIG. 6. (Color onhne) The optimal control pulses for (a) C{t) 
and (b) d{t) — do{t), where do{t) = 0.5t. The main features 
here are large perturbations at the initial and final time due to 
the boundaries, and slower modulations for the intermediate 
stage of the transport. 

the excitation at the initial time, and then 'decelerating' 
it near the final time. Small modulations are required 
at intermediate times in order to prevent the excitation 
from spreading over too many sites. It is also worth not- 
ing that the speed achieved here is at least two orders 
of magnitude faster than is possible in the adiabatic case 
for comparable fidelities [l[. 

If we decrease T, we find for all times T shorter than 
a particular time TqgL that even after applying the op- 
timisation algorithm we are still unable to achieve high- 
fidelity state transfer. In other words, there is a min- 
imum time required to perform the transfer (33| . The 
lower-bound on the value of Tqsl by the quantum 

speed limit (QSL); no transfer can take place faster than 
the QSL allows. Fig. [7] shows the same transfer of ex- 
citation as in Fig. [SJ but in this case we have set the 
total allowed time T = ^qg^ (how we determined Igg^ 
is shown later). One sees clearly that the evolution of 
the system is that of a wave of excitation, moving with 
an almost constant velocity along the chain. When we 
choose time T < Igg^, we find accordingly that the opti- 
mal control algorithm is unable find an optimal solution, 
even after many thousands of iterations. This is a strong 
indication that we have gone beyond the quantum speed 
limit, and there is no solution by which we can transfer 
the excitation across the chain in the given time. The 
evolution of the system in this case is shown in Fig [5] In 
comparison to Fig. [71 one sees that the evolution looks 
much the same. However, if one compares the excitation 
profile at T/2 for both evolutions, one sees that while the 
evolution at the QSL has the excitation wave centred at 
the 51st site (i.e. the halfway point), the evolution for 
a time T < Tq^^ falls short of the halfway point after 
T/2. This is an indication that we are indeed beyond the 
QSL, since if we cannot reach the halfway point before 
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FIG. 8. (Color online) The probability density of the wave- 
function along the chain at different times: (a) t = 0, (b) 
t = T/4, (c) t = T/2, (d) t = 3T/4, and (e) t = T, where 
T < Tqsl = 53.30J~^ Both d{t) and C{t) were found after 
100,000 iterations of the optimal control algorithm. 



half the time has elapsed, we might well guess that we 
cannot reach the final site in the remaining half of the 
time. 

We can see this failure of the optimisation algorithm 
more clearly in Fig. |9l For times T > IggL, the infi- 
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FIG. 9. (Color online) The decrease in infidelity of a transfer 
across a chain with 101 sites against the iterations of the 
control algorithm. The solid (red) line is the convergence for 
a transfer time T = 70.92J^^ > Tqsl, the dashed (green) 
line for a transfer time T = TqgL = 56.50J~^, and the dotted 
(blue) line for a transfer time of T < TqgL = 53.33J~^. 




FIG. 10. (Color online) The average speed of the excitation 
wave Vd versus Va- The solid (red) line shows the effect of 
filtering on the optimised pulses for a chain length of 41 sites, 
the long-dashed (green) line shows the same for 61 sites, the 
short-dashed (dark blue) line for 81 sites, the dotted (pink) 
line for 101 sites, and the dashed-dotted (light blue) line for 
121 sites. The black dotted line is the line Va = Vd- 



delity converges almost exponentially towards zero. For 
times T < I^qgL, the decrease in infidelity saturates after 
several hundred iterations. 

Another indication that the QSL has been reached can 
be found by examining the average "velocity" of the ex- 
citation wave as it moves across the chain. Given a 
total time T for the propagation, the average rate at 
which the excitation should be transmitted is given by 
Va = {N — i)/T. Examining the dynamics, we can see 
that for much of the propagation time, the excitation 
moves along the chain with an (approximately) constant 
velocity. We can quantify this velocity as 



Vd 



(12) 



where (x) = {'${t)\x\'i/{t)) is the expectation value of the 
position of the excitation along the chain. In other words, 
we take the average position of the excitation in the time 
interval [T/4, 3T/4] (to avoid effects at the ends of the 
chain) and divide by the average time taken to reach that 
position, T/2. 

In the ideal case, we would have Va ~ Vd, in which case 
the optimal solution would be the transit of the excita- 
tion along the chain at exactly the average rate required 
to reach the other end. However, as we cross the thresh- 
old set by the QSL, we should find that Vd reaches a 
maximum, which is the maximum speed at which the ex- 
citation can propagate. This is exactly what is seen in 
Fig. [101 

The last issue we want to address is robustness. In 
essence, how much information in the control pulses given 
in Fig. ini (and indeed in all of the control pulses at the 
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FIG. 11. (Color online) The solid (red) line shows the Fourier 
transform of the pulse d{t) — da{t) for a chain length of 101 
spins with a total time T — 56.50J^^. 




FIG. 12. (Color online) The infidelity of the transfer related 
to the maximum frequency component of the controls retained 
after filtering. For key, see Fig. (TO] 



QSL) can be discarded without detriment to the trans- 
fer fidelity? Figure [TT] shows an example spectrum of a 
pulse for d{t) for a transfer along a chain of 101 spins 
at the QSL, and Fig. [12] shows the effect on the fidelity 
after filtering the optimised pulses. The filter applied is 
a simple frequency cutoff: the pulse (in frequency space) 
is convoluted with a function 



liv; t'max) 



U if \u\ < I^max, 

otherwise. 



(13) 



where t'max is the maximum allowed frequency in the 
pulse. We see that not all of the frequencies in the con- 
trol pulses need be retained; on average, we only need 
frequencies up to around 4 J in order to maintain a high 
fidelity. Note that this is independent of the chain length 
N. Figure [T3| shows a set of pulses that transport the 
excitation along a chain of 101 spins with an infidelity 
/ < 10""*, where the maximum frequency component is 
- 4 J. 
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FIG. 13. (Color online) The optimal control pulses for (a) 
C(t) and (b) d{t)-do{t), where do{t) = 1.77*. The maximum 
frequency component is ~ 4J. 



IV. THEORETICAL LIMITS OF 
NON-RELATIVISTIC QUANTUM THEORY 

Does the limit Igg^ discussed in the last chapter have 
a physical origin, or is it simply a numerical constraint, 
stemming from the construction of the optimisation rou- 
tine itself? If, in fact, we are able to reach the physical 
limit by application of optimal control routines, then it 
would appear that optimal control can not only be used 
to improve the operation of experimental implementa- 
tions, but indeed to probe a system's dynamics and phys- 
ical limits. This connection was already investigated in 
Ref. [3^; here, we elucidate further the methods that 
were applied and more specific conclusions. 

The physical limits on quantum systems (and hence 
any physical system) have been investigated theoretically 
in quantum systems for several years; such considerations 
lead Lloyd in 2000 to calculate the maximum rate at 



which any machine can process information [32[ . In par- 
ticular, the notion of a "quantum speed limit" has been 
reported by several authors. We briefly recount this the- 
ory and its particular application to our problem. 



A. The quantum speed limit 

What is the absolute maximum speed at which we can 
transfer information along our chain? This amounts to 
finding the minimum time it takes for the given initial 
state |5'(0)) to evolve to the goal state \(Pn)- A possible 
route for finding this minimum time was explored by Car- 
lini et al. [33 |. where it was shown that one may derive 
the time-optimal Hamiltonian for a given state evolution 
by minimising the quantum action 5* of the system, by 
which the problem may be interpreted as a quantum ana- 



logue of the classical brachistochrone. In principle, the 
same procedure could be performed in our case, but the 
complexity of the calculation is prohibitive for a many- 
body system like ours. Hence we ask a somewhat simpler 
question, as in Giovannetti et al. [29| : how fast can a 
quantum system under a time-independent Hamiltonian 
evolve in time? 

The notions of energy and time are not inseparable, 
an idea that presents itself in the enigmatic time-energy 
uncertainty relation (45j . Hence the minimum time in 
which we can perform some given evolution must be con- 
nected to the related energy scales. This minimum time 
is referred to as the quantum speed limit (QSL). For the 
case where the evolution is from an initial state to an or- 
thogonal state for a time- independent Hamiltonian, this 
relation can be written explicitly as [soj 



'''QSL 



2AE J ' 



with 



AE ^ {m\[H{t) - Eimm) . 

E = {iiO)\H{tmO)) . 



(14) 



(15) 
(16) 



As pointed out, this is only valid when the time evolu- 
tion is governed by a time-independent Hamiltonian: E 
and AE are a measure of the energy resources available 
in the system only at the initial time, which for time- 
independent Hamiltonians defines a fixed energy scale. 
In our case, the methodology must be slightly modified, 
by considering instead the mean energy spread of our 
system as it evolves under our time-dependent Hamilto- 
nian, which we find by averaging the instantaneous en- 
ergy spread of the system over the time interval [0,r]. 
By integrating over time, we effectively apply the bound 
to infinitesimal time steps dt where the Hamiltonian 
is approximately constant. We modify the definition in 
Eq. dUl) to read pl 



TQSL = max < — 



2J' 2A£2 



where 



with 



A£^^j\E{t) dt, 



AEit)^^{^{t)\[H{t)-Em^'{t)), 
E{t) ^ {^{t)\H{t)\^5lit)) . 



(17) 



(18) 



(19) 
(20) 



As was already pointed out, this speed limit defines the 
time it takes to rotate from the initial state to an or- 
thogonal state. Since the initial and final sites are not 
directly coupled, we cannot immediately rotate from the 
initial state to our goal state. Due to this condition, we 
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postulate that the speed hmit must be interpreted as an 
effective time-per-site; the total time it takes to traverse 
the chain is this time-per-site multiplied by the number of 
sites (minus one) in the chain, or equivalently, the num- 
ber of edges we have between the initial state vertex and 
the final state vertex when one views the spin chain as a 
connected graph. 

Equation (|17p effectively states that the minimum time 
it takes to rotate from the current system state to an or- 
thogonal state is bounded from below by 7rfi./(2J) (we 
shall see later that for the evolutions we consider, the 
second term in Eq. ([T7)) is always less than this term, so 
that we can neglect it). By considering the speed limit 
of a simple two-spin system with a coupling strength J, 
we can associate this bound with the time it takes to 
swap an excitation between only two sites, given that for 
the initial state the excitation is completely localised on 
one of the two sites. Using the reasoning above, we see 
that the quantum speed limit theory predicts that the 
minimum time to traverse the chain is given simply by 
the time it takes to perform a swap between two neigh- 
bouring sites (which we shall henceforth refer to as 'or- 
thogonal swaps') mutiplied by the number of sites in the 
chain (minus one). However, in our particular system, 
at some intermediate time it may be (as we have already 
seen from the results in Section UlI)) that the excitation 
does not perform repeated swap operations, but rather 
moves along the chain as a dispersed "wave" . If one now 
imagines the picture of the excitation wave moving from 
site to site, we note that two excitation waves centred 
at neighbouring sites are not orthogonal, unlike when we 
have the excitation fully localised on a single site. This 
means that we can expect the actual propagation time 
to be shorter than the one calculated from simply doing 
repeated orthogonal swap operations. The optimised sys- 
tem performs a controlled excitation-wave propagation, 
which we can view as a cascade of effective swap oper- 
ations, each shorter in duration than that given by the 
orthogonal swap. We are then motivated to write the 
total time to traverse the chain as 



QSL 



7(iV - 1)tqsl , 



(21) 



where 7 is a dimensionless constant that quantifies the 
effective swap duration in terms of the orthogonal swap. 
As a side remark, we note that one can also imagine 
mapping the full chain with the effective swaps onto a 
shorter chain with orthogonal swaps, which is analogous 
to a reduction of the transmission length of the chain. 
Similar ideas have already been explored for long range 
interactions in Ref. |46|. 



B. Comparing limits 

As already alluded to in Section lTTll there comes a point 
where the optimal control algorithm is no longer able to 
reach an optimal solution. We aim to show that this 
limit on the evolution time (which we denoted by ^Qg^) 
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FIG. 14. (Colour online.) The infidelity reached after R = 
100, 000 iterations for different chain lengths A'' and eflfective 
time-per-site (T — b)/{N — 1), where b = 3.65. 



corresponds to the quantum speed limit for the system 
Tqsl discussed in Section HV Al 

The procedure for determining rgg^ is as follows. We 
select a chain length A'^, and set some initial evolution 
time T which we assume to be longer than the corre- 
sponding IgsL- We perform optimal control on the sys- 
tem for a fixed number of iterations R. We then repeat 
this for shorter and shorter times T. The results of the 
simulations are shown in Fig. [141 Note that we plot the 
effective time-per-site (T — b)/{N — 1) in order to make 
comparisons between chains of differing lengths easier. 
One sees clearly that for longer times, we arc able to 
complete the state transfers with high fidelities. As we 
reduce the time, we begin to see that the final value of 
the infidelity does not converge to zero, even after many 
thousands of iterations of the control algorithm. Some- 
where in between these two extremes lies the limit of the 
optimal control algorithm. We quantify this by setting a 
threshold e for the infidelity; the time ^qg^ for each N is 
defined as the smallest value of the time T for which the 
infidelity I < e after R iterations. This threshold obeys 
a linear relation; 

with a ~ 0.34 and 6 = 3.65. Note that this is a pos- 
teriori the same b used in the effective time-per-site 
(T - b)/{N - 1) for Fig. [H The introduction of the 
constant b describes additional effects due to the bound- 
aries of the chain, where the excitation wave is generated 
at the beginning of the evolution, and then collapsed into 
a localised excitation at the end. Additionally, b is not 
dependent on TV (unless the chain lenght is of the order 
of the width of the spin- wave) . 

We now compare the results from the quantum speed 
limit Tqsl with Igg^, which is shown in Fig. 1151 In order 
to evaluate Eq. (|17p for each value of N, we must numer- 
ically calculate the second term in the bracket, since it 
depends upon the time evolved state \ij{t)). For all points 
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FIG. 15. (Colour online.) A comparison of the quantum- 
speed-limit-time Tqsl with the optimal control limit Tqsl- 
The solid (red) line is Tqsl with 7 = 1, which is the repeated 
orthogonal swaps. The (blue) crosses are Tqsl for different 
A'^ in the ranse 21-131 with e = 5 ■ lO"'^. The dashed (green) 
line is Tqsl + b with 7 — 0.34. 

at our defined thresliold, this comes out to be less than 
the first term in the bracket in Eq. pT)) . so that the speed 
limit is given simply by the effective swap time. One finds 
that optimal control outperforms what can be achieved 
through applying repeated swap operations between ad- 
jacent spins. Furthermore, by ignoring boundary effects 
for T^qgL, we find that our model for the quantum speed 
limit fits the data with a value of 7 = 0.34. This means 
that the speed limit achieved with the optimal control 
can be described (ignoring the ends of the chain) as a 
cascade of effective swaps. 

V. CONCLUSION 

We have shown that we can successfully apply opti- 
mal control to the system given in Eq. ([2]) to produce 
fast transfers of excitations along spin chains; two or- 
ders of magnitude faster, in fact, than was reported in 
Rcf. [ij for comparable fidelities. This has application 
in the fast transport of quantum states over short dis- 
tances. Furthermore, we have found a fundamental limit 



for optimal control beyond which optimisation is not pos- 
sible, and identified it as a speed limit on the dynamics 
of the system, which is manifested by the dynamics as 
the propagation of an excitation wave with constant ve- 
locity. We compare this with the standard formulation 
of the quantum speed limit, and show that for our many- 
body problem, the quantum speed limit implies that the 
optimal strategy for transport is characterised by effec- 
tive swaps along the chain. We confirm this through a 
comparison with the numerical results. 

It is interesting to note that aside from the theory on 
the quantum speed limit, there is a large body of work 
concerned with a similar bound specific ally for spin sys- 
tems, namely the Lieb- Robinson bound j47H5(j |. It would 
be interesting to investigate the connection between this 
bound and the QSL, although it is likely difficult to quan- 
tify this explicitly. 

Wc have shown that not only is optimal control a use- 
ful tool for the optimisation of tasks relevant for quantum 
information processing (specifically transmission of quan- 
tum information along a spin chain) , but also as a means 
to probe the limits of many-body quantum systems where 
the theoretical methods become unwieldy. We expect 
that given the generality of the method, it should be 
able to probe fundamental limits of many quantum sys- 
tems that can be efficiently simulated. Indeed, wc used 
the same technique to prove a bound on the duration of 
a unitary SWAP operation on a spin chain, showing that 
it was achievable in a time that scaled only polynomially 
with the number of sites Q (although it was not shown 
that this was a fundamental limit). Wc will continue with 
such investigations in future work. 
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